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Abstract 

The recent progress in revealing classical integrable structures in quantum models solved by Bethe 
ansatz is reviewed. Fusion relations for eigenvalues of quantum transfer matrices can be written in the 
form of classical Hirota's bilinear difference equation. This equation is also known as the completely 
discretized version of the 2D Toda lattice. We explain how one obtains the specific quantum results 
by solving the classical equation. The auxiliary linear problem for the Hirota equation is shown to 
generalize Baxter's T-Q relation. 
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1 Introduction 



Acquiring some experience in classical and quantum integrable systems, one would probably say that 
these theories hardly have something to do with each other. They do look like different branches of 
mathematical physics with their own methods, notions and traditions. The correspondence is believed 
to exist but practically it does not seem to uncover itself to a satisfactory extent. 

Actually, the very questions used to be asked in the classical and quantum cases usually did not 
have common points. So there is no surprize that main ingredients of the theories practically do not 
intersect. For instance, the inverse scattering method has a very different meaning in the two theories 
to say nothing about such specific tools as the finite-gap integration technique (from the classical side) 
or the Bethe ansatz (from the quantum side). 

At the same time we would like to believe that such a strong property as integrability may add 
something important to the correspondence principle. At present it is gradually realized that there 
should exist a structure common for classical and quantum integrability, or rather that classical and 
quantum cases should be the two faces of this structure. 

Having all this in mind, let us give an outline of a more technical story presented here. These notes 
review the recent progress (see the first paper in ref. p]) in uncovering classical integrable structures in 
quantum models. More attention is paid to some questions which remained obscure or were not very 
clearly written in ffl . In particular, we try to follow, at least schematically, the whole long way from the 
simplest quantum i?-matrix to the classical discrete soliton equations. We hope that in the future the 
story may be considerably shortened and simplified. 

In Sect. 2 we review basic elements of the fusion procedure on the example of the simplest solution to 
the Yang-Baxter equation with rational dependence on the spectral parameter. The consisetency of the 
fusion procedure relies on the Yang-Baxter equation for the i?-matrix. As a result, one obtains a family 
of quantum transfer matrices (generating functions for commuting hamiltonians) acting in one and the 
same quantum space. They depend on some (spectral) parameters and commute for all values of these 
parameters. 

These transfer matrices are functionally dependent. They satisfy a number of functional relations 
(called fusion relations) which are analogues of relations for characters of linear groups but have a more 
complicated structure. Remarkably enough, this structure is well known in the classical soliton theory in 
a very different context. It is an integrable discretization of the 2D Toda lattice. In the bilinear form it 
is known as Hirota's bilinear difference equation (HBDE). This fact suggests an intriguing link between 
quantum and classical integrable systems lying much deeper than any kind of a naive " classical limit" . 
Let us stress that it is the quantum integrability (i.e. commutativity of the multiparameter family of 
quantum operators) that allows one to reduce operator relations to classical equations for eigenvalues. 
In this approach, different quantum states correspond to different solutions of the classical equation 
satisfying certain boundary and analytic conditions. 

Sect. 3 is a brief summary of main facts about HBDE. From the classical viewpoint, HBDE plays the 
role of a master equation for the majority of known continuous and discrete soliton equatuions. R.Hirota 
has shown that HBDE unifies various types of them: the Korteweg-de Vries (KdV) equation, Kadomtsev- 
Petviashvili (KP) equation, two-dimensional Toda lattice (2DTL), the sine-Gordon (SG) equation, etc 
as well as their discrete analogues can be obtained from it by different reductions, specifications of 
parameters and continuum limits. Furthermore, HBDE itself has been shown to possess soliton solu- 
tions and Backlund transformations. The key elements of the theory are discretized Zakharov-Shabat 
representation (the zero curvature condition) and the corresponding auxiliary linear problems. 

In Sect. 4, we show how the specific quantum notions and methods (such as the transfer matrix, 
Baxter's Q-operators, the nested Bethe ansatz, etc) are translated into the purely classical language. 
In other words, solutions of discrete classical nonlinear equations appear in the Bethe ansatz form, a 
surprizing fact that still waits for a deeper understanding. 
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2 Quantum fusion relations 

i?-matrix 

The fundamental rational i?-matrix acting in C k x C k has the form 

R{u)=u + 2V, (2.1) 

where V is the permutation operator, V{x®y) = y®x, and u is the spectral parameter. It is convenient 
to represent R{u) as k x k matrix in the first space (called the auxiliary space) with non-commutative 
entries which in their turn are k x k matrices acting in the second space C k (called the quantum space): 

k 

R(u) = ul + Eij ® Eji . (2.2) 

Here I is the unit matrix and is the k x k matrix such that (Eij) mn = Si m Sj n . 
This i?-matrix satisfies the Yang-Baxter equation 

i?i 2 (wi - U 2 )Ri3(ui ~ u 3 )R 23 (u 2 - u 3 ) = i?23("2 - u 3 )R 13 (ux - u 3 )Ri 2 {ui - u 2 ) . (2.3) 

Both sides are operators in C fc x C fe x C fe . We use the standard notation: R\ 2 {u) acts as R(u) in the 
tensor product of the first two spaces and as identity in the third one, similarly for R± 3 , R 23 . Graphically, 



eq. (2.3) can be written as follows: 




which means that the line u 3 can be moved over the intersection point of the other two lines. 
T-matrix 

The quantum monodromy matrix (T-matrix) is the following product of i?-matrices in the auxiliary 
space Vo = C k : 

T(u) = R QN (u - y N ) . . .R 02 {u - y 2 )Roi(u - yi) . (2.4) 

This is a k x k matrix with operator entries. They act in the quantum space <8>£Li Vi, V\ = C k . Matrix 
elements of i?oi and Roj for i ^ j commute with each other because they operate in different components 
of the tensor product. The parameters yi are arbitrary; sometimes they are called rapidities. 



The following basic relation is a direct corollary of eq. (2.3): 

R\ 2 {u - v)T 13 (u)T 23 (v) = T 23 (v)T 13 (u)R 12 (u - v) , (2.5) 

where the two auxiliary spaces are Vi = V 2 = C k and the quantum space is V3 = ®f = iVi. 
Graphically the T-matrix looks as follows: 
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Traditionally, the auxiliary space is attached to the horizontal line while the quantum space is assigned 
to the vertical lines with rapidities yi. 



Fusion procedure 



Using the i?-matrix (2.1) as a building block, it is possible to construct more complicated solutions to the 
Yang-Baxter equation. This is done by multiplying the fundamental solutions and subsequent projection 
on irreducible representations of GL{k). We give the general scheme of the fusion in the auxiliary space 
which is of prime importance for our purposes. 



Let A = (Ai, A2, . 
YlTLi = n - Let 



, A m ) be a Young diagram with n boxes and m lines Ai > A2 > 

-Pa : ®} =l Vi -> 



> A m , m<k, 



be the projection operator on the space of the irreducible representation of GL(k) corresponding to A. 
Write in the box with coordinates (i-th line and j-th column) the number 



S(ij-) =u-2(i- j) . 



(2.6) 



For example: 



u 


u + 2 


u + 4 


u + 6 


11 — 2 


u 


u + 2 




u — 4 







(2.7) 

Enumerating all the boxes in the natural order (from left to right in the first line, then continue from 
left to right in the second line and so on), one gets the ordered sequence of numbers s\, S2, . . . , s n (in 
the example above s\ — u, S2 = u + 2, S3 = u + 4, S4 = u + 6, S5 = u — 2, sq = u, s-j = u + 2, sg = u — 4). 

The i?-matrix acting in <g> C fc is 



R (x \u) = P x R nQ (s n ) 



R2o(s 2 ) ®Rio(si)P\ 



(2.8) 



The fundamental i?-matrices Rio(si) are tensor multiplied in the auxiliary spaces Vi = C , their entries 
being multiplied in the common quantum space Vq = C k in the indicated order. Graphically, 



Px 





Sn 












Si 







Pa 
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The important property of R^ (u) is given by 
Theorem 2.1 R^(u) satisfies the Yang-Baxter equation in Vi ® V2 ® V3: 

- R 12 ) ( M l _ "2)i?i3 ) ( w l - "3)^23(^2 ~ U3) = ^23(^2 - lia)-?^ 0"l - U3)i?^''(Ml - u 2 ) , (2.9) 

w/iere Vi = V" (A \ V 2 = V 3 = C fc . 

The idea of the proof is simple: to represent the projector P\ as a product of fundamental i?-matrices. 



This representation alows one to move the lines using the Yang-Baxter equation (2.3) for fundamental 
i?-matrices. However, implementation of this idea in general case is quite involved from technical point 
of view. We find it useful to outline main steps of the proof since it simultaneously helps to recover the 
structure of zeros of the fused i?-matrix, which will be important later. 



The following two elementary properties of the fundamental R- matrix (2.1) are crucial for the con- 
struction: 

a) There exists a value of u at which R(u) is proportional to the permutation operator: 

flia(Q) = 2Pi 2 . (2.10) 

Graphicaly this means that if spectral parameters atached to a pair of crossed lines coincide, one can 
"eliminate the intersection point" as follows: 



U 

#12(0) = 2 - 



u 



b) The i?-matrix degenerates at the points ±2: 

J? 12 (±2)=4P 1 ± 2 , (2.11) 
where P^ progects on symmetric (antisymmetric) tensors in Vl®V%. 

Recall the ordered sequence of spectral parameters 81,3%, . . . , s n introduced at the beginning of this 
section. The property b) suggests to try 

n 

P ~^Y[ RlJ ( Sl - Sj ), (2.12) 

i<j 

where the structure of the ordered product corresponds to the pattern of complete intersection of n 
lines carrying spectral parameters Si . The order of successive intersections is irrelevant by virtue of the 
Yang-Baxter equation. For example (n — 3): 
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One can prove by induction that this prescription does work for "hook diagrams" A = (Ai, 1, 1, . . . , 1) 
(i.e. such that any diagonal contains exactly one box) providing the desired projection operator. 

However, in general case we encounter a difficulty: p™ alve can be identically zero. Indeed, if A — (2, 2) 
we have 




that is zero since the complementary projectors P + and P~ meet together. A little inspection shows 
that this zero arises for the same reason each time when there are at least two lines i ^ j such that 
Si = Sj. Such lines always exist if the Young diagram contains a diagonal with at least two boxes. In 
this case a sort of regularization is necessary. 

Let us introduce a small parameter e and modify the sequence si, s%, . . . , s n writing 

s\ ij) =u-2{i-j) + {i-l)e (2.13) 
in the box with the same linear order s{ , s|, . . . , s e n . Consider the operator 

n 

^ e) =n^-( s i- s i) ( 2 - 14 ) 



with the same order in the product as in eq. (2.12). 
Lemma 2.1 As e — > 0, 

p| e) = C e K P x + 0(e K+1 ) , (2.15) 

where n = #{(ii, Ji), (*2 5 J2) |*i < *2> il ~ ji — *2 — J2} is the number of pairs of boxes on the diagonal 
lines and c is a numerical constant. 

This lemma is sufficient to prove the Yang-Baxter equation for R^ (it) . 

The fusion procedure in the quantum space is made in a similar way. This procedure allows one 
to define the i?-matrix iJWW(u) acting in ® for two arbitrary Young diagrams A, fi. The 
Yang-Baxter equation then reads 

- u 2 )R^(u 1 - u 3 )R { ^(u 2 - U3 ) = R^(u 2 - u 3 )Rg\ Ul - u a )B$ M (u, - u 2 ) , (2.16) 
where Vi = V< A \ V 2 = V 3 = C fc . 
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Each matrix element of (u) is a polynomial in u. However, the fusion procedure brings to R^ (u) 
a number of "trivial" zeros in the sense that they are common for all the matrix elements. In other 
words, there are some values of u at which R^ (u) is the zero operator. These zeros should be extracted. 

Using the same argument as above, it is easy to see that 

R w {u- Si ) = 0, i = 2,3,...,n, (2.17) 

the remaining n-th zero being non-trivial. If Sj = Sj for i ^ j, we have zero of higher degree. Extracting 
the common multiplier, we redefine 

RW(u)^R^(u)(j[s^ (2.18) 

that is a polynomial of degree 1 with operator coefficients. 

At last, we specialize this formula for rectangular Young diagrams of length s and height a: A = 
(s, s, . . . , s) (a times): 

(o-ls-l \ _1 /a-l \ _1 

lH[(u-2q + 2p)\ (U(u-2q)\ . (2.19) 
g=0p=l ) \q=l J 

For a = 1 the last product should be skiped, for s = 1 the double product should be skipped. 



Fused T-matrices 



Similarly to eq. (2.4) one can introduce the T-matrix with the auxiliary space Vo — and the quantum 
space 



T^(u) = R^(u - y N ) . ..R$(u - y 2 )Rhi ) (u - Vl ) . 



(2.20) 
(2.21) 



It follows from fl2.16[ ) and fl2.20| ) that 

fl8 )M (« - f^'M^W = T^{v)T£\u)R^\u v) , 

where the two auxiliary spaces are Vi = V^ x \ V2 = V**', and the quantum space is V3 = ®f— 1 Vi. 

Each factor in ( |2.20[ ) gives an independent contribution to the common polynomial multiplier. Let 



us explicitly extract this multiplier in the case of rectangular Young diagrams a x s: 

/ a— 1 s—1 



T (axs) (m) = T (axs) ( u ) J] JJ 0( u - 2 q + 2 P ) 

\g=0p=l / 



a-l \ 

n ^ u - 2 «) ■ 



Here 



N 



<t>(u) = Y[(u-Vj) 



(2.22) 



(2.23) 



and r T^ axs '(u) denotes the renormalized T-matrix. For generic polynomial <fi(u) it contains non-trivial 
"operator" zeros only. If there are pairs i,j such that y^ — yj = 2, T^ aXs '{u) acquires some extra 
"trivial" zeros and should be renormalized further. This case corresponds to higher representations in 
the quantum space. 



Quantum transfer matrices 

The key notion of the theory is quantum transfer matrid^ obtained by taking trace of T^(u) in the 
auxiliary space: 

T (x \u) = Tr aux T^ ( u _ A x + K) , (2.24) 

1 The name "transfer matrix" came from statistical mechanics on the lattice. It is a bit misleading since the transfer 
matrix is a scalar in the auxiliary space. The word "matrix" here is related to the quantum space. 
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where A' denotes the transposed diagram (i.e. reflected with respect to the main diagonal), so X[ is 
height of the first column of A. The shift of the spectral parameter is introduced for later convenience. 

The crucial property of quantum transfer matrices is their commutativity for all values of the con- 
tinuous spectral parameter u and the discrete spectral parameter A: 

[rW(tt),TW(tj)] = 0, (2.25) 



This fact immediately follows from the Yang-Baxter relation (2.21) 



Therefore, we have a family of commuting operators in the quantum space. The point is that 
hamiltonians of quantum integrable systems belong to this family, so T^(u) are integrals of motion. In 
other words, (it) is a generating function for commuting hamiltonians. Simultaneous diagonalization 
of transfer matrices T^'(u) is one of the corner-stones of the theory. In this way one obtains spectral 
characteristics of quantum integrable systems. 

For the case of rectangular diagrams a x s we introduce the special notation 

T s a (u) = T^ axs \u) (2.26) 

which will be used in the sequel. 



Extracting the "trivial" zeros (2.22) amounts to changing the normalization: 



(n-ls-l a-1 \ 1 

Y[ ]]_(f>{u-s-a + 2l + 2p + 2) Y[(/)(u-s-a + 2l)) . (2.27) 
1=0 p=l 1=1 ) 

Each eigenvalue of the transfer matrix is a polynomial in u and all the " trivial" zeros (common for all 
the eigenvalues) are removed. 

Quantum determinant 

The maximal "external power" of the fundamental T-matrix T{u) in the auxiliary space C k is called 
the quantum determinant of T(u). More precisely, we put 

D {k) {u) = detT(u) ee THu) = T {lk) {u -l + k). (2.28) 
i 

The last equality is due to the fact that the fused auxiliary space is one-dimensional. In other words, 
D^ k \u) is a scalar in the auxiliary space by definition. The following theorem shows that D^ k \u) is a 
scalar in the quantum space too. 

Theorem 2.2 The quantum determinant D^ k \u) of the fundamental T-matrix T(u) with k- dimensional 
auxiliary space lies in the center of the algebra generated by matrix elements ofT{u): 

[£><*>(«), {T(v)) ij ] = 0, i.j 1.2 /,-. (2.29) 

It is given explicitly by the formula 

fe-i 

£)( fc ) ( u ) = <f>(u + k + 1) JJ 4>(u - k - 1 + 21) , (2.30) 

i=i 



where 4>{u) is defined in eq. (2.25). 



The idea of the proof is to use the Yang-Baxter equation which allows one to reduce everything to the 
quantum determinant of the fundamental i?-matrix. In this case the assertion can be verified by a direct 
computation. 



Let us show how to verify eq. (2.30) without any computation. The last k — 1 factors contain just 



the "trivial" zeros coming from eq. (2.27), so we only need to explain the first factor. Note that at u = 2 



the fused i?-matrix contains the following block: 
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s k = -2k- 

s 2 = -4 
81 = -2 




which is zero since it is the projector on the ( fc + l )-th external power of C fc . So i?< lfc )(2) = as the 
whole thing. This yields the first factor in eq. ( [2.30 ). 



There is a more general formula derived by the same argument, 

/fe-ls-l \ fc-l 



21 + 2p + 2) JJ 4>(u - s - k + 21) , T°(u) = 1 . (2.31) 



\l=0 p=l 

which will be used in Sect. 4. 
Functional relations 



A combination of the fusion procedure and the Yang-Baxter equation results in numerous functional 
relations (fusion rules) for quantum transfer matrices. 

We illustrate the origin of these relations on the simplest example. Consider the product 



Tl(u - 2)T 1 1 (u) = Try l8 y 2 (T 20 (u - 2) ® T 10 (u)) , 



(2.32) 



where Vi = V2 = C k are auxiliary spaces and the quantum space Vq is arbitrary. Insert the unit matrix 
/ = P + + P~ represented as sum of the complementary projectors inside the trace. Using cyclicity of 
the trace, the property (P^) 2 = P^ and the Yang-Baxter equation, we arrive at the relation 



Tl{u - 2)T£(u) = Tf{u - 1) + T 2 \u - 1) . 
A more general relation derived in a similar way is 

T> + l)Tl(u -8) = T^=^(u) + Ti +1 (u) 



(2.33) 
(2.34) 



Proceeding further, it is possible to show that T^ x '(u) can be expressed through either T^(u) or 
Tf(u) only. 



Theorem 2.3 The following determinant formulas hold true: 

TW(u) = det {Tl_ l+J {u - Ai + Ai - A* + - 1)) 
IS 2 jS^i 



where it is implied that TW(tt) = 1 (% denotes the empty diagram). 



TW(«) = det (Tt- i+j (u X[ + Ai + \> 

1<ij<Ai 



(2.35) 
(2.36) 



We remind th e read er that A' de notes the transposed diagram. The entries T^u) (resp., T[ a (u)) of the 
matrix in eq. (2.35) (resp., (2.36)) are transfer matrices corresponding to the one-line (resp., one-column) 
diagrams. 



9 



These formulas were obtained by V.Bazhanov and N.Rcshctikhin. Sometimes they are called quantum 
Jacobi-Trudi formulas. They are " Yang-Baxterization" of the classical Jacobi-Trudi identities in the sense 
that the latter are reproduced from the former by forgetting the dependence on u. 



The specification of eqs. ( |2.35 ), (2.36) to rectangular diagrams reads 

T«(u)=det(T} +l _ J (u-i-j+a + l)), i,j = l,...,a, !?(«) = 1, (2.37) 

T s » = det (T? +i - j (u -i-j + 8 + 1)), i,j = l,...,s, T » = 1 . (2.38) 

These determinant formulas exhibit a very nice structure. Though, they give only a partial solution 
because T^(u) or Tf(u) (entering as "initial data") are still to be determined. 



Bilinear form of the fusion rules 



The fusion rules can be recast into another suggestive form which is in a sense more practically useful 
since it allows one to obtain a complete solution. At the same time, this form provides a link with 
classical non-linear integrable equations, which is our main concern here. 

It follows from eqs. ( 2.37 ), ( 2.38| ) that transfer matrices for rectangular Young diagrams obey a closed 
set of relations among themselves. Using the Jacobi identity for determinants, they can be represented 
in the following model-independent bilinear form: 



T:(u + 1)T S > - 1) - T^T^u) = T: +1 (u)Tr\u) 



(2.39) 



Since T"(u) commute at different u, a, s„ the same equation holds for all eigenvalues of the transfer 
matrix, so from now on we can (and will) treat as a number- valued function. Note that the 3 

variables u, s, a enter almost symmetricallyin spite of their very different nature. 



Let us mention an analogue of (2.39) for more complicated diagrams. Consider, for example, Young 
diagrams consisting of two rectangular blocks (i.e. with a\ lines of length s\ + S2 and the rest <Z2 lines 
of length si) and let T^f^(u) be the corresponding transfer matrix. Then it holds 



t: 



a\ — l,a2 



Sl + l,fl2 



t: 



0-1,0,2 — 1 



Sl+l,fl2 



(u+1)T ;™ (u _ 1) 

(u + i)t;^+> - 1) + t:i;^(u + i)2£+-i (u i) . 



(2.40) 



In what follows we mainly deal with eq. (2.39) only. 



Equation ( 2.39 ) is the top of the quantum construction described in this section. Remarkably, this 
very equation is well known in the theory of classical soliton equations. This is famous Hirota's bilinear 
difference equation (HBDE) for the function T"(u) of 3 variables. The reason of this coincidence is 
hidden somewhere in the representation theory. We believe that there should exist a formulation of 
the quantum theory such that this coincidence is obvious from the very beginning. However, we are 
not going to dwell upon this aspect of the problem. Our aim here is more technical. In the sequel, we 
treat the fundamental relation ( [2.39] ) not as an identity but as a fundamental equation and show how 
to extract from it the specific quantum information. The solution to (classical!) HBDE then appears in 
the form of Bethe equations. In our opinion, this fact is quite remarkable by itself. On the other hand, 
we anticipate that this approach makes it possible to use some specific tools of classical integarbility in 
quantum problems. 



Comments and references 

1. Throughout this section we use notations and terminology of the quantum inverse scattering method 
(QISM) developed by the former Leningrad school. At present QISM (= algebraic Bethe ansatz) provides 
the most convenient framework to analyse the phenomenon of quantum integrability. The early paper 
remains one of the best reviews of the subject; see also the book ||. The graphical interpretation is 
based on the reformulation in terms of the factorized scattetring . 
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The fusion procedure was invented by P.Kulish, N.Reshetikhin and E.Sklyanin g, Q and generalized 
to elliptic solutions of the Yang-Baxter equation by I.Cherednik 0. The regularization ( [2.13 ) was 
suggested in ref. ]^], for the general proof of the analogue of Lemma 1.1 in the elliptic case see |J. 

The notion of quantum determinant was introduced in the paper ]l6| for 2x2 monodromy matrices 
and generalized to k x k matrices in the paper j^]. For the complete proof of Theorem 1.2 in the elliptic 
case see e.g. [jL7j. 

2. Functional relations for transfer matrices have a long story. Some of them appeared for the first 
time in the papers ||, |10fl . The general quantum Jacobi-Trudi formulas ( 2.35| ), (2.36) appeared in Q , 
see also JT^ ]. The bilinear form of functional relations was first suggested by A.Klumper and P.Pearce in 
the particular case of Ai-type models and then generalized to the Afc_i-case by A.Kuniba, T.Nakanishi 
and J.Suzuki fl4j|. They also suggested bilinear fusion relations for models associated to other Dynkin 
graphs. Eq. ( 2.40| ) an d its ge neral izations to other types of Yang diagrams can be obtained from 
determinant formulas ( [2.35| ), ( 2.36 ) by means of the Pliicker relations. 

3. The transfer matrices play the role of quantum characters. They may be also called "quantum 
Schur functions" . Emphasizing this analogy, we use notations and conventions from the book jl5| on 
symmetric functions. 

4. The results described in this section are valid for elliptic sol utio ns of the Yang-Baxter equation, 
too (sometimes with minor changes). Instead of Yang's i?-matrix (2J2) one should start from Belavin's 
i?-matrix |L8] with elliptic dependence on the spectral parameter: 



where 



R{u)=Y,K^ u ) E "'® E i 

(?]( U + 2)\kT) 



3-1 
k 



k ' 2 



(2r)\kr)G 



if i + j = i' + j'(mod k) and R^j, (u) = otherwise. Here 



3-1 
k 



(2.41) 



(2.42) 



(i]u fcr) 



a 






(u\ 



lez 



exp 



(iiTT(l + a) 2 + 2iir(l + a)(u + 0)) 



(2.43) 



is the Jacobi theta-function with characteristics. There are two parameters: 77 and r. (This form of 
Belavin's i?-matrix is taken from refs. ]l9| ], p0[.) 

Again, the i?-matrix degenerates at u = ±2: R(—2) = M_P~, R{+2) = P + M + , where M± are some 
invertible matrices. The caref ul ana lysis shows that the fusion procedure goes through leading to the 
same basic functional relation ( 2.39 ). Formula ( 2.31 ) for the quantum determinant also holds with the 
following function </>(u): 

1 



N 



(u) = Ip 



(2.44) 



We call functions of this form elliptic polynomials in u. The number TV is degree of the elliptic polynomial. 



The elliptic solution ( 2.41 ) is no longer GL(fc)-invariant but is still associated with the root system 
Ak-\. The monodromy matrices built out of this i?-matrix form representations of generalized Sklyanin 
algebras, 2-parameter deformations of U{glk) [§l), p2|- We refer to integrable models with the i?-matrix 
( [2.41 ) as Ak-i-type models. 



3 A classical view on Hirota's equation 

This section is a brief survey of HBDE in the purely classical context. Firstly we list various types 
of HBDE and show their equivalence. Then discrete versions of the zero curvature representation and 
auxiliary linear problems are presented and Backlund transformations are discussed. 
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At the first glance, the content of this section has nothing to do with specific tools of the algebraic 
Bcthe ansatz partially described in the previous section. Indeed, the two sections may be read indepen- 
dently of each other. However, in the next section we link them together and show that main elements 
of the quantum theory can be translated into the language of classical integrability and embedded into 
the purely classical context. 

Equivalent forms of the bilinear equation 

A) Hirota 's original form. In Hirota's original notation it is 

(zi exppi) + z 2 exp(D 2 ) + z 3 exp{D 3 j) T ■ T = , (3.1) 

where Zi are arbitrary constants, T — T(x±, x 2 , x 3 ) is a function of 3 variables and Hirota's D-operator 
Di = D Xi is defined by 

F(D x )f(x)-g(x)=F(d v )f(x + y)g(x-y) . (3.2) 

y=o 



In the more explicit notation eq. (3.1) looks as follows: 

ziT(x x + l)T(xi - 1) + z 2 T(x 2 + l)T(x 2 - 1) + z 3 T(x 3 + l)T(x 3 - 1) = (3.3) 

(here and below we often skip variables that do not undergo shifts). 

Note that the 3 variables enter in a symmetric fashion and the equation is invariant under their 
permutations (and a simultaneous permutation of Zj's). The equation is also invariant under changing 
the sign of any one of the variables and under the transformation 

T(xx,x 2 , x 3 ) -> Xq(%1 +x 2 + x 3 )xi(x 2 +x 3 - xi)\2(xi + x 3 - x 2 )x 3 {xi + x 2 - x 3 )T(xx,x 2 ,x 3 ) , (3.4) 

where Xi are arbitrary functions. 
The transformation 

T(x u x 2 ,x 3 ) -> z~ x?/2 (-z 2 )-^/ 2 z 3 ^ /2 r(x 1 , x 2l x 3 ) (3.5) 



converts eq. (3.3) into the canonical form, 

T{x x + l)T{xi - I) - T(x 2 + l)T(x 2 - 1) + T(x 3 + l)T(x 3 - 1) = (3.6) 

which does not contain any free parameters. 

After appropriate identification of the variables one recognizes the bilinear fusion relation ( 2.3£ ) . We 
refer to x±,x 2l x 3 as direct variables (in contrast to light cone ones, see below). 

A') "Gauge invariant" form: 

v( . , sw -. \ (l + y(ari,a;2,a:3 + l))(l + y(a;i,ar2,a;3-l)) ,„ 

Y(x u x 2 + l,x 3 )Y(x u x 2 - l,x 3 ) = (1 + y - 1( , i + M2i , 3))(1 + y - 1(xi _ 1; ^, 3)r (3-7) 

1 (Xl + 1,X2,X 3 )1 (Xl - l,x 2 ,x 3 ) 



where 



is a gauge invariant quantity: the "gauge" transformation (3.4) does not change it 



B) KP-like form. The equation for a function r(pi,p2,p 3 ) of 3 variables reads 

Zi t( Pi + l) T (pa + 1,P3 + 1) + z 2 t{ P2 + l)r(pi + l,p 3 + 1) + z 3 r{p 3 + 1)t( Pi + l,p 2 + 1) = . (3.9) 
Again, one can eliminate the arbitrary constants by the transformation 

t(pi,P2,Ps) f - — J (~^J T bl'-P2,P3) ■ (3.10) 



12 



C) 2DTL-like form: 

r x (t, t+ l)T B (t + 1, t) - r x (t, t)T s (t + 1, i + 1) - rr x+1 (t, t+ l)r a _i(t + 1, t) , (3.11) 

where T x (i, f) is a function of the 3 variables and r is an arbitrary constant. The variables t, i are called 
light cone coordinates. Note that in this form the permutation symmetry is lost. However, an analogue 
of eq. (3.4) holds true: 

r x (t, t) -> X o(2x + 2t)xi(2t) X 2(2t)X3(2x - 2t)T x (t,t) 
sends solutions to solutions. The transformation 

T x {t,t) ^r- x *l\ x {t,t) 

sends eq. ( fo.ll ) to its canonical form (r = 1). 

At last we present the linear substitutions which make canonical forms of equations a), b), c) equiv- 
alent. 



(3.12) 



(3.13) 



A)<-» B): 



B)<-> C): 



A)«-» C): 



T(xi,X 2 ,X 3 ) = T(pi,P2,P3), 
Xl = P2 + P3 , X 2 = Pl + p 3 , 2^3 = Pi + P2 ; 

t(pi,P2,P3) = T x (t,t), 

Pi = t, p 2 = x-i, P3 = t; 

T(Xi,X2,X 3 ) = T x (t,t), 

x± = x , x 2 =x + t~i, x^ — t + i; 



(3.14) 



(3.15) 



(3.16) 



Clearly, these linear substitutions are not unique. All other possibilities can be obtained from the 
given one by applying a transfo rmat io n of t he form (x\, X2, x$) — * (±Xpn), ±ip(2), ±Xp(x\), where P is 
a permutation. Using formulas ( 3.14 )-( 3.16 ) one can easily obtain gauge invariant forms of equations B) 
and C). 



Discrete Zakharov-Shabat representation 

The reformulation of classical nonlinear integrable equations as flatness conditions for a two-dimensional 
connection is the basic ingredient of the theory. The flatness means that subsequent shifts along any pair 
of the time flows commute. These conditions are known as Zakharov-Shabat equations or zero curvature 
representation. 

Let us consider a family of difference operators acting in the space of scalar functions of a variable u: 

MW = -z+ 1 ; f j e- a " , (3.18) 

where z, z are arbitrary parameters, the shift operator e du obeys e ±9 "/(w) = f(u ± l)e ±a ", r(u) is 
a function of u depending also on a number of extra variables t,p, . . . ;i,p, . . . When it is necessary to 
show the dependence on these variables explicitly, we write, e.g. MW(i,p, . . .), etc. We call M {z \ AfW 
M -operators; note that their coefficients are gauge invariant. 
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Theorem 3.1 The discrete Zakharov-Shabat equations (zero curvature conditions) 

MW(p,t + i)M&(p,t) = M (z) (j) + l,i)M w (|),(), 

AfW(p,i+ l)M^(p,i) = M<*)(p + l,t)M< ffl >(p,t) , 
M (2) (t + 1, t)M {z) (t, t) = M (z) (t, t+ 1)P (t, t) 
are equivalent to the following bilinear relations for t{u) = j- > < p, P'"'(u): 



r P+l,ti 



zt- ' ^(u)r p ' t+1 (u + 1) - wt p <* +1 (uK +M (« + 1) + # x (p, t)T p+ ^ t+1 {u)T p - t {u + 1) = , 
zr p+1 '*(u + 1)t p ' I+1 (u) - wt p ' I+1 (u + 1)t p+1 '*(u) + H 2 (p, t) T p+1 ' i+1 (u + 1)t p '\u) = , 
zT^ +1 (u)T t+u (u) + H 3 (t,t)T^(uy+^ +1 (u) = (^)-V' t + 1 (u + l)r t + 1 ' t "(u-l), 
respectively, with arbitrary u-independent functions Hi. 

The proof consists in straightforward commutation of the M-operators. 
We see that 



(3.19) 

(3.20) 
(3.21) 

(3.22) 

(3.23) 
(3.24) 



and 



eq. ( |3.22|) coincides with HBDE in the KP-like form (3J3) under the identification T p,t (u) = r(p, t, u) 



z = zi, 



w 



-Z2, 



Hi= z 3 ; 



and 



eq. ( |3.23 ) coincides with HBDE in the KP-like form (3J) under the identification r p '*(u) = r(p, t, —u) 



w 



-Z2, 



Ho 



Z3 ; 



and 



eq. fl3.24|) coincides with HBDE in the 2DTL-like form (|3.1lj ) under the identification t*'*(u) = r tt (t, t) 



The variables t,i provide a family of commuting flows parametrized by z, z. In fact one can realize 
M-operators as difference operators in any one of them. The variable in which M-operators act will be 
called the reference variable. The functions Hi are automatically fixed to be the constants if one requires 
the simultaneous compatibility of all the corresponding zero curvature conditions. 



Linearization of HBDE 

The zero curvature conditions are equivalent to compatibility of an overdetermined system of linear 
difference equations. These linear equations are called auxiliary linear problems (ALP). They play a 
very important role in the theory. Common solutions to ALP (called wave functions) carry the complete 
information about solutions to the nonlinear equations. All the properties of the latter can be translated 
into the language of the ALP. This is what we mean by the linearization of HBDE. 

The discrete Zakharov-Shabat equations ( 3.19 )-( ^21 ) imply the compatibility of the linear problems 



M (z V <*(u) = -znj) t+l ' l {u) , (3.25) 

M^^(u) = -z^ i+1 {u) (3.26) 

for any discrete flows labeled by z, z. Note that the "eigenvalues" in the r.h.s. can be arbitrary: this is 
the matter of redefinition of ?/> ' («). Our choice corresponds to the smooth continuum limit. 

More explicitly, eqs. (Uf), ( p^ ) read (see ftTTj\ ), (pj|)): 

zV*'*(u)^*'*(u) - + 1) = z^ t+1J {u) , (3.27) 

z^' 1 («) - C*'* (u)^*' f (u - 1) = I/' * +1 («) , (3-28) 

where 



t'''(ii)t w ''(«|1) 

T*+M"(u) T t,t( u+ l)' 



(3.29) 
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C*'*(u) 



l (u-iy- i+1 (u + i) 



r t,t+i 



(u)t*'*(u) 



These formulas become more symmetric in terms of the " urmormalized" wave function 

p(u) = ip(u)r{u) . 



Making this substitution in (3.25), (3.28), we get: 



ZT 



t+M 



{u + l)p^{u) 



{U)p tr \u + 1) = ZT^{ U + l) P t+1 '\u) , 



r t,t+l 



r t,t+l 



(u + l)p £ '> - 1) = zV^Wp'-'+V) 



(3.30) 

(3.31) 

(3.32) 
(3.33) 



Duality 



The ALP ( 3.32| ), (3.33) have a remarkable property: they are symmetric under interchanging r and 
p. Furthermore, one may treat them as linear problems for the function r, the compatibility condition 
being a nonlinear equation for p. This equation is again HBDE. We refer to this fact as the duality 
between "potentials" r and "wave functions" p. 



More precisely, rewriting eqs. ( 3.27 ), ( 3.2S ) as linear equations for 

„*+i,t+i( u+1 ) 



ir («) = 



t+M+i 



(u + 1) 



- 'V> t+M+1 



(u + 1) 



we have 



-<)„ 



z)V^{u) ) ^\u) = -^ t " 1 '*(u) , 



'z - C*- f (« + l)e 9 " ) = # t,f_1 («) , 



(3.34) 
(3.35) 



where V" and (5 are given by the same formulas (3.29), ( |3.30 ) with p in place of r. After the replacing 
p — » r the difference operators in the l.h.s. become formally adjoint to the operators (3.17), ( 3.18| ) 
(defined by the rule (f(u)e kdu Y — e~ kdu f(u)). It then follows that the compatibility conditions arc 
described by Theorem 3.1 with r replaced by p. 

Therefore, passing from a given solution r to p we have got a new solution to HBDE. This is a 
Backlund-type transformation. 



Backlund flows 

One may repeat the procedure described above once again starting from p and, moreover, consider a 
chain of successive transformations of this kind. Let us introduce an additional variable m to mark steps 
of the "flow" along this chain and let r^'(tt) , Pm{u) be r and p at m-th step. 

The Backlund flow is defined by 

p^) = A^W. (3.36) 

This mean s tha t r at t he next step of the " Backlund time" m is put equal to a solution p of the linear 
equations (3.32), (3.33) (up to the factor z u ). Then these linear problems become bilinear equations for 



Jk+l 



(uy m+1 (u + 1) - T £ l (u + i)tI +1 { u ) + + ly+iM = , 



where t (resp., i) in eq. (||7|) (resp., ( p.38| )) is skipped. 

Similarly, defining the second Backlund flow (the Backlund time is now denoted by fh), 



Pm{ u ) = (z) u t^ +1 (u+ 1) 



we get from ( p2] ), ( |3~33| ): 



r t+i 



M^+iW - iWi+il") = ( zz ) 



.ifu + ljr^u-l), 



(3.37) 
(3.38) 



(3.39) 



(3.40) 
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rUuy^iu + 1) - + 1) + r^{u + lH l+1 (u) = . (3.41) 

In these equations one immediately recognizes different forms of HBDE described by Theorem 3.1. So, 
we conclude that the Backlund flows can be identified with the commuting discrete flows. 

Comments and references 

1. The bilinear difference equation was suggested by R.Hirota in the paper p3| which summarized his 



earlier studies on discretization of nonlinear integrable equations J24j]- j28j. HBDE (3.6) can be viewed 
as an integrable discrete analogue of the 2-dimensional Toda lattice. 

2. As it was first noticed by T.Miwa J29|, discrete Hirota's equations can be obtained from the 
continuous KP hierarchy by choosing the time flows to be certain infinite combinations of the standard 
continuous flows of the hierarchy (Miwa's transformation). This approach was further developed in the 
papers ^] , |3l| . The different forms of HBDE arise when one applies Miwa's transformation to different 
continuous hierarchies. The function T u (t,t) can be identified with the r- function of the hierarchies 
(restricted to a finite number of the discrete flows). 



3. The "gauge invariant" form (3.7) of HBDE is a discrete counterpart of nonlinear integrable 
equations written in terms of potentials and fields rather than r-functions. Eq. (|3.8| ) is a discrete version 
of the famous formula U(x) = 2&i logr(a;). Some particular cases of equation (|3.7|) emerge naturally in 
thermodynamic Bethe ansatz |32[-[p4|. 



4. An example of the discretized zero curvature representation for eq. (3.6) was given by R.Hirota 
p3|. In the physical language, the discrete connection is a lattice gauge field. The approach emphasizing 
the relation to gauge field theories on the lattice was developed by S.Saito and N.Saitoh |56|. We have 
presented these results in a modified form which makes the theory completely parallel to the 2DTL 
theory. The bilinear form of Backlund transformations was discussed by R.Hirota ^?J. The r <-> p 
duality was pointed out in rcf . p6[ . 

5. Similarly to the 2DTL hierarchy, there exists an infinite hierarchy of bilinear difference equations. 
For the explicit form of higher equations of the hierarchy see ref. pqj. We remark that the general 



quantum fusion rules in the bilinear form (see e.g. eq. (2.40)) coincide with the higher HBDE-likc 
equations (after a linear change of variables). 

6. In the continuum limit {z — > oo, z — > oc) eqs. (3.27), ( |3.28 ) turn into the familiar ALP for the 
2DTL f37|: 

Otlpn = tpn+1 + d t {\0g ^^-)lpn , 

d^n = T " +1 a"~V n-i , (3.42) 

T 

where n is identified with u. The compatibility condition is 

dtT„diT n - T n d t dtT n = t„ + it„_i , (3.43) 

which in terms of 

lp n (t,t) =l0g — -=- 

T n (t,t) 

acquires the form 

Olefin = e^~^-i - e^+^» , (3.44) 
the first non-trivial equation of the 2DTL. 

4 Bethe ansatz results from the Hirota equation 

We are going to show how to reproduce the quantum Bethe ansatz results by solving the classical Hirota 
difference equation ( 2.39| ): 



t> + i)r s > - 1) - r^^dw = Tr»ir>) ■ (4.i) 
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In general, HBDE has many solutions of very different nature. To extract the solutions of interest, we 
should first of all specify boundary and analytic conditions for T s a (u). 



Boundary conditions in a and s 

The values of T s °(u) and Tjf(u) for A^-i-type models should be considered as boundar y con ditions (b.c). 
Recall that T°(u) = 1 and T^(u) is the quantum determinant given explicitly by eq. ( 2.31 ): 

s-i / k—i \ 

Ts(u) = [| U{u - s + k + 2 + 2p) Yl 4>{u - s - k + 2p + 21) J . (4.2) 
p=o \ 1 = 1 ) 

Let us note that T^(u) obeys the discrete d'Alembert equation: 

T*(u + l)T*(u - 1) = T* +1 (u)Tl x (u). (4.3) 



This is easily seen from eq. (4.2). To put it differently, recall that the general solution to the discrete 
d'Alembert equation is x+( u + s )X- ( u ~ s )> i- e - it is factorized into a product of a " holomorphic" function 
X+ and an " antiholomorphic" function \-- It is possible to represent T*(u) in this form indeed. For 
instance, in the rational case we have: 

Xk{u - s) 

Xk (u) = 2^ ft (r(^pi + 1) ]]r(^ + o) , (4.4) 

i=l \ 1=1 J 

where T(u) is the gamma-function. Combining this property with the Hirota equation, we arrive at the 
following b.c: 

T s °(it) = as a < or a > k. (4.5) 

The b.c. in s reads 

T s a (u) = as - k < s < 0, and < a < k . (4.6) 
The meaning of this condition will be discussed later. 



Analytic conditions in u 



A very important condition in models on finite lattices (which follows, eventually, from the Yang-Baxter 
equat ion) is that T^{u) for any fixed a, s, u has to be an elliptic polynomial in the spectral parameter u 
(see ( 2.44 ) multiplied possibly by an exponential function. 

More precisely, the gauge invariant quantity 



Y s a (u) 



^i(")r,"-i(t») 

T s a+1 («)T s a - 1 ( U ) 



(4.7) 



has to be an elliptic function of u having 27V zeros and 2N poles in the fundamental domain. This 
implies that T®(u) has the general form 



T s » = A a s e 



N(a,s) 

n < 



1/2 
1/2 



(4.8) 



where 



„(«.») 



, Ag, /i(a, s) do not depend on u (and meet some additional requirements). For models with 



rational R- matrices T® (u) degenerates to a usual polynomial (multiplied by the exponential function). 
We note that the b.c. (4.5), (|4.6|) are consistent with the required analytic conditions. 
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Normalizations 



By normalization we mean multiplying T"(u) by the "gauge" factors Xoi u + s + a )Xi( u + s — a )X2( M — 
s + a)x3(w — s — a). Boundary and analytic conditions may depend on the gauge. We point out the 
following three distinguished normalizations. 

i) Det-normalization was already discussed in detail: T®{u) — 1, T^(u) is the quantum determinant 



ii) Minimal polynomial (MP) normalization. The" gauge" invariancc (3.4) allows one to remove all 



zeros from the characteristics a± s ±it = const. These are just the "trivial" zeros brought by the fusion 



procedure. The minimal polynomial appears in the gauge (2.27). The boundary values at a = 0, k then 
become: 

T s » = ftu-a), 

T s fe (u) = <f>(u + s + k). (4.9) 

Note that in this normalization T° is an " antiholomorphic" function (i.e. depends on u + s only) while 
T h is a "holomorphic" function (i.e. depends on u — s only). The most important feature of the MP 
normalization is that the elliptic polynomials T"(u) have one and the same degree N for all values a, s 
(when Tg(u) is not identically zero). In the sequel, we use the MP normalization. 

iii) Canonical normalization. Since the functions T s °(u) and Tjf(u) at the "boundaries" obey the 
discrete d'Alembert equation, they can be gauged away, that is can be made equal to unity: T®(u) = 
Tj?(u) = 1. This allows one to simplify the equations for the price of imposing much more complicated 
analytic properties of solutions. Here we will not discuss this normalization. 

Zero curvature representation 

Let us represent the bilinear fusion relation (f4.l|) as a discrete zero curvature condition. To do that, we 



can make use of Theorem 2.1. Consider eq. (3.24) and identify 

!?(«)= r*- f (tO, (4.10) 

where 

s = -t-t, a = u + t-i (4.11) 
(cf. ( |3.16| )). We have to substitute d t — > — d s + d a , <% — > — d s — d a and 



t,t=const 



a.s- 



=const 



9 a , 



so the M-operators read 



r°( M )r°+ 2 ( M + i) 

'T s °+/(u)T s 0+1 (u+l)' 



'•) = «^ a °- ,, - (4.12) 

(4.13) 



M a.s = -^H 1 e " ' 

1 ' ' T« = l{u)T S {u) 



In terms of the operators 



M = e d °- da M(a, s), M = e d ° +da M(a, s) (4.14) 



the Zakharov-Shabat equation (3.21) acquires the form of a commutativity condition: 

[M, M]=0. (4.15) 
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Linear problems 



The commutativity of two operators implies the existence of a common eigenfunction: 

M^ as (u) = E^ as (u), M^ as {u) = E^ a ' s (u). 
Explicitly, these equations read 



#°' s (u + l) -z- 



T s U(u)T?(u + l) 



Passing to the " unnormalized" wave function, 

F?(u)=T«(u)y a ' s (u), 

we get 

T s U(u)F s a (u + 1) - zT^{u + l)Fr\u) = ET?{u + lj^u) , 
T^u + l)Fr\u - 1) - zT:iI{u)F:{u) = ET" (u) F"Zi (u) . 



(4.16) 

(4.17) 
(4.18) 

(4.19) 



(4.20) 



Now, to identify eq. ( 3.24 ) with eq. (LI) literally, we set z ~ z ~ = —1. Besides, redefining 
^ a ' s {u) — > E( s ~ a ^ 2 (E)( s+a ^ 2l >i> a ' s (u), we can always choose E = E = 1 without loss of generality. In 
this way we get the following ALP: 

T££(u)F?(u) T?+\u + l)F« +1 (u 1) = T-^F^iu) , 

T« +1 (u + l)F s a (u) - T«{u)F« +l {u + 1) = T«+\u + l)F£*(u) . (4.21) 

An advantage of the light cone coordinates t, t is that they are separated in the linear problems 
(compare ( p2j ), ( [3~33| ) with (ggg)). However, in contrast to a, s,u they do not have an immediate 
physical meaning. 



Due to the duality property (Sect. 3) F^(u) obeys the same HBDE: 

Ff(u + l)F«(u - 1) - F^F^iu) = F? +1 (u)F?~ l (u) . 



(4.22) 



We require F^(u) to have the same analytic properties as T^(u) (though, degree of the elliptic polynomial 
may be different). 



The b.c. (L5) allows one to impose a similar condition for F: 

F s a (u) — as a < or a > k — 1 



(4.23) 



so that the number of non-zero functions F is one less than the number of T's. The functions F a at the 
ends of the Dynkin graph (a = 0, k — 1) have a very special form. From the second equation of the pair 
(4.21) at a — and from the first one at a = k — 1 it follows that F®(u) (respectively, F*~ depends 
on one light cone variable u — s (resp., u + s). We introduce a special notation for them: 



F»(u) = Q k - X {u - s), F*-\u) = Q k - X (u + s) 



(4.24) 



Furthermore, it can be shown that the important condition (L9) relates the functions Q and Q as follows: 

Qfc_i(«)=Q fc _i(u + fe-l). (4.25) 

Therefore, the analytic properties and b.c. for F®(u) are the same as for T^{u) under a substitution 
4>{u) by Qk-i{u). The only change is a reduction of the Dynkin graph: k — » k — 1. Using this property, 
one can successively reduce the -problem up to A\. Below we use this trick to derive Ak-i ("nested") 
Bethe ansatz equations. 
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Nested Bethe ansatz and Backlund flows 



Quantum integrable models with internal degrees of freedom can be solved by the nested (hierarchical) 
Bethe ansatz method. The method consists essentially in integration over a part of degrees of freedom 
by an ansatz of Bethe type, the effective hamiltonian being again integrable. Repeating this step several 
times, one reduces the model to an integrable model without internal degrees of freedom which is solved 
by the usual Bethe ansatz. 

The classical face of this scheme is a chain of Backlund transformations, i.e. passing from solutions 
of the non-linear equation to (properly normalized) solutions of the auxiliary linear problems discussed 
in Sect. 3. 

To elaborate the chain of these transformations, let m = 0, 1, . . . , k mark steps of the flow Au-i — > A± 
and let F s a,m (it) be a solution to the linear problem at (fc — m)-th level. In this notation, F£' k (u) = T^(u) 
and F s a ' fe ~ 1 (u) = F£(u) is the corresponding wave function. For each level m the function F^' m (u) obeys 
HBDE of the form (|4 22j) with the b.c. 



F^ m (u) =0 as a < or a > m . 



(4.26) 



The first (a = 0) and the last (a = k — 1) components of the vector F^' m (u) obey the discrete d'Alembert 
equation and under the condition (4.24) are " antiholomorphic" and " holomorphic" functions respectively. 
We denote them as follows: 



(u) = Q m (u - s) . 



l (u) = Q m (u + s), 



where it is implied that Qk(u) = 4>(u). Furthermore, it can be shown that the relation 

Q m (u) = Q m (u + m) 
can be imposed simultaneously for all 1 < m < k. 



The linear problems (4.21) at level m, 



f: 



a+l,m+l 



s+1 



(u)F^ m (u)-F l 



a+l,m+l / 



(u + l)F^(u - 1) = F^(u)F^ m (u) , 



f:^ +i (u+i)f:^(u)-f: 



(4.27) 
(4.28) 



(4.29) 
(4.30) 



l (u)F^(u + 1) = + l)F£? m (u) 

look as bilinear equations for a functions of 4 variables. However, eq. (4.29) (resp., eq. ( f.30| )) leaves 
the hyperplane u + s — a = const (resp., u — s — a = const) invariant, and actually depends on three 
variables. 



Restricting the variables in eq. (4.29) to the hyperplane u + s — a = 0, i.e. setting 



T u {t,a) = F° 



1 ( \ 

(u) 



(4.31) 



we reduce it to the 2DTL-like form of HBDE (3.11) for T u (m,a) with m and a being the light cone 
coordinates. The b.c. is 



T„(m, 0) = Q k -m(2u), T u (m, k - m) = Q k _ m (m - k) = const. 



(4.32) 



A similar equation can be obtained from the second linear problem (4.30). They are nothing else than 
eqs. ( |3.38D , fl3.4l| ) for the Backlund flows. 

It is convenient to visualize this array of r-functions on a diagram; here is an example for the A3-case 
(fc = 4): 






1 
















Qx(u - 


>) 


Qi(u + s) 












Q2{u - 


>) 


F^(u) 


Q 2 (w + s) 










Qz{u - 


s) 


F^(u) 




Q 3 (u + s) 








Qa{u - 


s) 


T s( u ) 




T?(u) 


Qa{u + s) 



(4.33) 
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Functions in each horizontal (constant m) slice satisfy HBDE, whereas functions on the u — s — a = const 
slice satisfy HBDE with m, a being the light cone variables. 

In the nested Bethe ansatz approach, the functions F£' m (u) are auxiliary objects - eigenvalues of 
transfer matrices for "intermediate" models arising at (k — m)-th step (level) of the hierarchical Bcthc 
ansatz^. The functions Q m {u) play a distinguished role. They can be identified with generalized Baxter's 
Q-operators in the diagonal representation (see below). In the general solution to HBDE these functions 
are arbitrary functional parameters. The additional requirement of ellipticity determines them through 
the Bethe equations for their zeros. 



Bethe equations as a discrete dynamical system 



Recall that the function r u (m,a) = F^L U m ( u ) (4-31) obeys HBDE in light cone variables: 

T u (m + 1, a)T u (m 7 a + 1) - t„(to, a)T u (m + 1, a + 1) = t u _i(to + 1, a)T u+1 (m, a + 1) . (4.34) 

Since t u (to, 0) = Qfc_ m (2u), nested Bethe ansatz equations can be understood as "equations of motions" 
for zeros of Q m (u) in discrete time m (level of the Bcthc ansatz). The simplest way to derive them is to 
consider the auxiliary linear problems for eq. (4.34). Here we present an example of this derivation in 
the simplest possible form. 

Let us assume that Q m (u) has the form 



M„ 



Q m (u) = e^ u [] 8 



1/2 
1/2 



(V(u-u?)) 



(4.35) 



(note that we allow the number of roots M m to depend on m) . Since we are interested in dynamics in m 
at a fixed a, it is sufficient to consider only the first linear equation ( 3.32) ) substituting r*'*(u) — » T u (m, a). 
An elementary way to derive equations of motion for roots of r«(m, 0) is to put u equal to the roots of 
p m,0 (u), p m ' n (u + 1) and p m+1 '°(u), so that only two terms in ( |3.32| ) would survive. Combining relations 
obtained in this way, one can eliminate p's and obtain the system of equations 



Q m -i{u™ + 2)Q m (uf - 2)Q m+1 (u] 1 ) 

Qm-l«)Qm(«7* + 2)Q m +l{< - 2) 



(4.36) 



as the necessary conditions for solutions of the form ( 4.35 ) to exist. With the "boundary conditions" 



Qo(u) = 1, 



)k(u) = 4>(u), 



(4.37) 



this system of Mi + M2 + . . . + Mk-i equations is equivalent to the nested Bethe ansatz equations 
for Ak-i-type quantum integrable models with Belavin's elliptic i?-matrix. The same equations can be 
obtained for the right edge of the diagram ( 4.33j ) from the second linear equation. In what follows we 
explicitly identify our Q's with similar objects known from the Bethe ansatz solution. 



Let us remark that the origin of equations (4.36) suggests to call them as equations of motion for the 
elliptic Ruijscnaars-Schneider (RS) model in discrete time. Taking the continuum limit in m (provided 
M m = M does not depend on to), one can check that eqs. (4.36) do yield the equations of motion for 
the elliptic RS model with M particles. The additional limiting procedure 77 — > with finite r/Uj = Xj 
yields the well known equations of motion for the elliptic Calogero-Moser system of particles. 

However, integrable systems of particles in discrete time have a richer structure than their continuous 
counterparts. In particular, the total number of particles may depend on the discrete time. Such a 
phenomenon is possible in continuous time models only for singular solutions, when particles can move 
to infinity or merge to an other within a finite period of time. Remarkably, this appears to be the case for 
the solutions to eq. ( 4.36 ) corresponding to eigenstates of quantum models. It is known that the number 
of excitations at TO-th level of the nested Bethe ansatz solution does depend on to. In other words, the 
number of "particles" in the associated discrete time RS system is not conserved. At the same time the 

1 are eigenvalues of 



2 Unlike eq. (4.1), eqs. (1.29), (1.30) can not be understood in the operator sense since F 
operators acting in different quantum spaces. 



and F n 
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numbers M m may not be arbitrary. It can be shown that for models with elliptic i?-matrices in case 
of general position M m — (N/k)m, where N is the number of sites of the lattice (degree of the elliptic 
polynomial <fi(u)). In trigonometric and rational cases the conditions on M m become less restrictive but 
still these numbers may not be equal to each other. 

Difference equations for Q m 's 

The functions Q m {u) obey certain linear difference equations. In principal, they can be obtained from 



the system of linear problems (4.29), (4.30) at levels m = 1, 2, ... k by excluding all functions F except 
those on the boundaries of the array (4.33). For k = 2 this can be done without any problem. Indeed, 
in this case there are only two non-trivial linear equations: 

Tj +1 (u)Qi(u - s) - T}{u + l)Qi(« -s-2) = 4>(u- s)Qi{u + 8+1), 

Tl +l {u + l)Qi(u + s) - T 1 s {u)Q 1 {u + s + 2) = cf>{u + s + 3)Q 1 (u - s - 1). (4.38) 

Dividing both sides of the first equation by <f>{u — s) and making use of the fact that the r.h.s. does not 
depend on u — s, one arrives at 

<t>{u - s)Q 1 {u - s + 2) + (f>(u - s + 2)Qx{u - s - 2) = A(u)Qi(u - s), (4.39) 

where 

T}(u+1) 

Recall that the b.c. is Tl 1 (u) = 0. Whence one obtains famous Baxter's relation 

0(u)Qi(« + 2) + 4>{u + 2)Q x (u - 2) = Tl{u)Q x {u) (4.40) 

which is a 2-nd order difference equation for Q\. The second linear equation yields the same result 
(recall that Qi(u) = Qx{u + 1)). 

In general the procedure becomes quite involved. Here we only quote the results: 

k 

J2(-l) a T?(u - a + l)Qx[u - 2a + 2) = , (4.41) 

a=0 

A Tf( M + + 1) Qk-i(u + 2a) 

^ ' <j>{u + 2a + 2) <f>(u + 2a) ' [ ' ' 

The equations for Q m 's with 2 < m < k — 2 have a more complicated form. 

Factorization formulas 

At last, we are to identify our Q m 's with Q m 's from the usual nested Bethe ansatz solution. This is 



achieved by factorization of the difference operators in (4.41) and (4.42) in terms of Q m {u). One can 



prove the following factorization formulas (looking like discrete Miura transformation): 

k 



-2d u Qk(u)Qk-l(u + 2)\ ( e -2d u _ Q2(u)Ql(u + 2) \ f-2d n _ Ql( U ) ) (443) 



Q k (u + 2)Q k - 1 (u) J V Q 2 (u + 2)Q 1 (u) J \ Q^u + 2 

V^/ ]) a-k T l( u + a + 1 ) 2ad u 

^ 1 4>{u + 2a + 2) 

a=0 v ' 

29„ Qiiu) \ ( 2du Q2(u)Qi(u + 2)\ ( 2du Q k (u)Q k - 1 (u + 2)" 



Qx{u + 2))\ Q 2 (u + 2)Q 1 {u)J "\ Q k (u + 2)Q k -i{u)J ' l ' ; 
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Note that these operators are adjoint to each other. The factors in the r.h.s. resemble M-operators 



Let us note that eq. (4.44) can be rewritten as follows: 



* Tf(u + g-l) 2at 

^ y ' 0(u + 2a-2) 



a=0 



Qi(u + 2) 28t 

1 Qi(u) 



1 _ Q 2 {u + 2)Q 1 (u-2) ^ 

Q2{u)Ql{u) 



1 _Q fe (n + 2)Q fc - 1 (,-2) e29 A 
Qk(u)Qk-i(u) ) 



Now the factors in the r.h.s. resemble M-operators ( 3.18| ). It is possible to show that coefficients of the 
operator inverse to (4.45) give T\{u): 



E 



T*(u + s-l) 2sdu 
<f>{u) 

Q k (u + 2)Q k _ 1 (u-2) 2du 



i{u + 2)Q 1 (u-2) 2du 



I _ e 

}k(u)Q k -i(u) J V Q2(")Qi(u) 7 V 



Qi(u + 2) 



2a,. 



(4.46) 



These formulas yield T-^u), T}(u) in terms of elliptic polynomials Q m with roots constrained by the 
nested Bethe ansatz equations which ensure cancellation of poles in T\ a (w). The transfer matrices T"(u) 
for a, s > 1 can be then found with the help of determinant formulas (2.37), @.3£ ). 



Comments and references 



1. With the b.c. (4.5) HBDE (4.1) is known as the bilinear form of the discrete two-dimensional 
Toda molecule equation |38| (the Toda lattice with open boundaries), an integrable discretization of the 
conformal Toda field theory [^9|. In particular, at k = 2 we have a discrete analogue of the Liouville 
equation. 



The b.c. in s of the form (4.6) is known in the classical theory, too. The b.c. of this kind is a 
hall-mark of forced ("semi-infinite") hierarchies of non- linear integrable equations emerging naturally in 
matrix models of 2D gravity (see e.g. J40[]). For instance, in the forced 2DTL hierarchy, the r-function 
T n (ti,t 2 , ■ ■ ■ | i i , 1 2 5 ■ ■ ■) is equal to zero at n — — 1: r_i = for any 



Similarly, imposing the analytic condition on the T-function of the form (4.8) is a familiar story in 
classical nonlinear integrable equations since the paper J4"i| , where the elliptic solutions to the KdV 
equation were studied. A systematic approach to elliptic solutions of the KP equation based on the 
finite-gap integration methods was developed by I.Krichever fl42| . 

So we see that each one of the boundary and analytic conditions has been known in the classical 
theory. However, they never met altogether. The specifics of solutions relevant to quantum problems is 
perhaps just in the combination of the above conditions which looks quite unusual for the soliton theory. 

2. The RS model (in continuous time) was introduced in the paper [|43"[ as a relativistic generalization 
of the Calogero-Moser system of particles. This model was shown to be integrable and the Lax repre- 
sentation was found. Recently, it was shown 1 45 that the RS system describes the dynamics of poles of 
elliptic solutions (zeros of r-function) to the 2DTL. Equations of motions for the discrete time analogue 
of the RS model were written down in ref. |Q, where an ansatz for the Lax pair was suggested. The 
close connection with the nested Bethe ansatz was observed and explained in ref. [0 . 

3. The fact that in the elliptic case degree of the elliptic polynomial Q m (u) is equal to M m — (N/k)m 
(provided rj is incommensurable with the lattice spanned by the two complex periods l,r and N is 
divisible by k) follows directly from Bethe equations (4.36). Indeed, the elliptic polynomial form of 



Q m (u) implies that if u™ is a zero of Q m (u), i.e., Q m (u™) — 0, then u™ + ni + n 2 T for all integers 
nx,ri2 are its zeros too. Taking into account the well known monodromy properties of the ^-function, 
one concludes that this is possible if and only if 



A I 



m+l 



2M„ 



(4.47) 
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which has a unique solution M m — (N/k)m satisfying the required b.c. This means that the nested 
Bethe ansatz scheme for elliptic Afc_i-type models is consistent only if N is divisible by k. 

In trigonometric and rational cases the conditions on degrees of Q m 's become less restrictive since 
some of the roots can be located at infinity. The equality in the formula for M m becomes an inequality: 
M m < (N/k)m. A more detailed analysis shows that the following inequalities also hold: 2M\ < M2, 
2M 2 < Mi + M 3 , ■ . ., 2M m < M m _i + M m+1 , N = M k > 2M k ^ - M fe _ 2 . 

4. It should be noted that the family of commuting transfer matrices generally does not define a 
quantum system uniquely. To do that, one should choose a hamiltonian, i.e. take a particular represen- 
tative of the commuting family of operators. So, to introduce a quantum integrable model to be solved, 
one has to give a prescription how to get hamiltonian from the transfer matrix. To be definite, let as 
consider the Ai-case. The hamiltonians are taken to be linear in logT s , to wit 

if = £ a p log T s (4.48) 

p 

where u^ p ' are some particular values of the spectral parameter. (For example, the hamiltonian of the 
spin-1/2 Heisenberg magnet is H = d u logTi(u) at u = 0.) If the hamiltonian is local, one can define the 
elementary energy function e(u) end represent the total energy as E — £(%•); where the momenta Uj 
of " quasiparticlcs" satisfy Bethe equations. This very important characteristics of the quantum system 
depends on the choice of H . 

Is there a room for e(u) on the classical side? The answer seems to be in the affirmative. This function 
might encode a way in which potentials or fields entering soliton equations are expressed through the 
r-function of a given hierarchy. (In the most popular example of KdV this is U(x) = 29^ logr(x).) 



Since T s (u) has been identified with a r-function, eq. (4.48) resembles such expressions in the hierarchies 
of discrete soliton equations. The hamoltonians H might then be identified with the corersponding 
potentials (or rather logarithms of potentials in the discrete case). 



5. The l.h.s. of eq. ( 4.45| ) is known as the generating function for T*(u). These formulas for the 



generating function coincide with the ones known in the literature (see e.g. |llj, |47|, [|48|]). They 
may be considered as non-commutative analogues of generating functions for symmetric multivariablc 
polynomials (see Remark 3 to Sect. 2). 
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